Local constitutive parameters of metamaterials 
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The reason of the non-locality of constitutive (material) parameters extracted in a usual way from 
the reflection-transmission coefficients of composite slab at moderately low frequencies is explained. 
The physical meaning of these parameters is clarified. Local constitutive parameters of metamaterial 
lattices are discussed and their existence at moderate frequencies is demonstrated. It is shown how 
to extract local material parameters from the dispersion characteristics of an infinite lattice and 
from reflection and transmission coefficients of metamaterial layers. 
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I. INTRODUCTION 

In the modern scientific literature the so-called 
metamaterials^ correspond to a lot of exciting results 
and a lot of mess. Metamaterials (MTM) are often 
defined as media which possess constitutive (material) 
parameters not observed in nature. Thanks to these 
material properties, many exotic phenomena were pre- 
dicted and discovered in the MTM structures (a detailed 
overview is presented ir>i). As a rule, MTM are presented 
by lattices of reciprocal optically small resonant particles, 
such as complex-shape metal inclusions (in different fre- 
quency ranges), plasmonic and polaritonic nano-spheres 
and nano- wires (at infrared waves and in the visible). 
Particles characteristic size 6 and maximal period of lat- 
tices a at the frequency of particle resonance are usually 
assumed to be much smaller than the wavelength A in the 
lattice matrix. This assumption is usually considered as 
allowing one to homogenize MTM and to explain their ex- 
otic properties through material parameters (MP). How- 
ever, it is important that i5 and a in MTM though small 
are not negligible with respect to A. Practically, they lie 
within the frequency band 



0.01 < 



A 



< 0.1. 



(1) 



As a result, the effective wavelength in MTM (it can be 
expressed though the lattice eigenwave wavenumber q as 
Acff = 27r/g) being shortened at the particle resonance 
compared to A can approach to a/ 2. At such frequencies 
the lattice spatial resonance (often called as Bragg 's res- 
onance) holds. Then the Bragg (staggered) mode which 
has complex wavenumber q = n/a -\- jlm(q) can be ex- 
cited in the lattice. Within the band of the staggered 
mode the homogenization of the metamaterial lattice is 
meaningless. If one formally introduces material parame- 
ters for such a regime they cannot satisfy to the physical 
conditions of locality. Basic physical limitations which 
are equivalent to the concept of locality are as follows 
(see e.g. irj^): 

• passivity (for the temporal dependence e~*"* it im- 
plies Im(£) > and Im(/i) > simultaneously at 



all frequencies, for e-''^* the sign of both Im(£) and 
Im(/i) should be negative); 

• causality (in the region of negligible losses the fre- 
quency derivatives of both Re(e) and Re(/x) should 
be positive); 

• absence of radiation losses (it was postulated for 
periodic infinite structures in^, then derived in^, 
now this principle is often attributed to workSi). 

Constitutive parameters of MTM are non-local at certain 
frequencies inside the resonant band of the particle. From 
it one often deduces that it is impossible to introduce 
local MP over the whole resonant band. In the present 
paper we will show that this opinion is wrong. 

It is evident that the staggered mode of a lattice does 
not cover the whole region of moderately low frequencies 
n and there are bands where a < Acff /2. There we should 
expect the locality of measured or simulated MP if they 
are measured or simulated properly). However, inspect- 
ing well-known works^&SiiiiiSiiiii'^ (and many others) 
devoted to the extraction of material parameters (MP) 
from measured or simulated reflection (R) and transmis- 
sion (T) coefficients of a metamaterial slab one observes 
another result. At least one of two extracted MP in all 
these works violates all locality conditions over the whole 
frequency rangeQ] It will be shown below that the reason 
of it is the different physical meaning of these MP than 
the meaning of local constitutive parameters introduced 
in the quasistatic theory of lattices (e.g. — ). Ini& it was 
properly noticed that MP introduced for orthorhombic 
lattices iniLiS are non-local by definition. In fact, these 
MP are exactly the same as MP measured or simulated 
jjj7,8j9jiiji2ji3ji4? ^ jj-^ -j-j^g present paper these non-local 

MP will be complemented by local ones extracted from 
same R and T for the same frequency range defined by 

m 

Ini one proposed to extract MP of lattices of thick- 
ness d as if these were filled by a continuous medium 
with parameters £off, ^efi- Then R and T for the normal 
incidence take form: 



R 



r{l - e-2j9cttd) 
1 — r^e^2-'9"tf'' ' 



„ j.2p-2jq^!id 



(2) 
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Here QcS = w^eoA*oecffMcff = koy/ScSfJ-cS is the wavenum- 
ber of the medium fiUing the layer. Also in these formulas 
r is the reflection coefficient from semi-infinite medium 
expressed through Z (medium wave impedance normal- 
ized to that of free space) as 



Z+1' 



z = 



II. 



LOCAL AND NON-LOCAL MATERIAL 
PARAMERETS 



In works22i22iS4i2& it was shown that the approximation 
of a uniform slab filled by effective continuous medium 
is too rough within the frequency range Q This is so 
even for non-resonant inclusions studied and, 
moreover, it is so for resonant ones. However, this obser- 
vation does not mean that the homogenization of slabs 
at moderately low frequencies is impossible at all. In 
these works two slightly different algorithms of homog- 
enization were developed. The original composite slab 
of thickness d was replaced by a three-layer structure. 
These effective layers are assumed to be filled by differ- 
ent continuous media. The central layer was described 
through local constitutive parameters El and /xi, which 
were calculated through the known polarizabilities of par- 
ticles. Two transition layers of small thickness dt ~ 5 
were described by special anisotropic material parame- 
ters Et and /it (also local). For Et and /if the closed form 
relations were derived. These parameters allow one to 
match normal components of averaged electric displace- 
ment < D > and magnetic flux < B > at all boundaries. 
This turned out to be possible up to qa ~ 1 and for 
all angles of incidence^. The last fact means that MP 
of both central and transition layers at fixed frequency 
do not depend on the direction of propagation. It is an 
obvious feature of local MP. Using this model the the- 
oretical refiection coefficient for arbitrary number N of 
lattice unit cells across the layer (starting from = 1) 
turned out to be practically equal to that obtained by 
exact simulationsSiS^. 

Though the use of local MPi22i2i2iiS^ makes the prob- 
lem of the plane- wave reflection in the frequency range 
more difficult these parameters are, at least, as important 
as non-local MP. Once extracted, tensors of local MP are 
independent on the wave incidence angle and polariza- 
tion. Non-local MP at the same frequency should be mea- 
sured or simulated separately for all possible directions of 
the wave vector. Local MP can describe the interaction 
between the medium with wave packages, moreover, with 
evanescent waves. Non-local MP are applicable for the 
only propagating wave. For the only wave they do not 
violate the passivity condition in the effective medium. 
In a plane wave electric and magnetic fields are related 
through the wave impedance. Then the negative elec- 
tric or magnetic losses (wrong sign of EcS for magnetic 
MTMl-ii4 or, vice versa, wrong sign of /Xeff for dielec- 
tric MTM) are compensated by the positive magnetic or 



electric losses. A similar speculation can be done on the 
causality. For wave packages and for evanescent waves it 
is not so because in this case the electric field and mag- 
netic field can be locally decoupled. For example, in the 
near zone of a wire antenna its electric field dominates. 
Then the wrong sign of electric losses in the substrate 
implies the energy generation, and the effective medium 
turns out to be active. It is also possible to show the vi- 
olation of the causality in these cases. So, non-local MP 
give not enough insight on MTM, and to know local MP 
is important. 

In the present paper it will be shown how to find lo- 
cal constitutive parameters through R and T of a slab. 
However, before it, let us understand why MP obtained 
i n Z — and other similar works are non-local. To answer 
this question one can discuss the meaning of non-local 
MP in terms of the transmission-line theory. Consider 
a problem of the plane-wave reflection by a semi-infinite 
lattice of dipole scatterers. The strict analytical solu- 
tion of this problem was found in^. For simplicity let us 
restrict by the special case of normal incidence to a or- 
thorhombic lattice of dipoles located in free space. Also 
we assume that the frequency range satisfies to the con- 
dition fcpa = u>ay/EofJ-o < TT, where a is the lattice period 
in the direction of propagation (in this frequency region 
the only eigenwave can propagate). Then the reflection 
coefficient of the lattice takes form^^S.: 

(fco-g)a 

(4) 



sm ■ 



sm 



{ko+q)a 



H. 



Here q is the eigenwave wavenumber and the factor H for 
which a closed-form expression was obtained irjSI takes 
into account all polaritons excited at the interface. It was 
proved in^*' that at the frequencies of our interest this 
factor has unit absolute value. So, it can be represented 
as an imaginary exponential H = e-'*. Then the influence 
of polaritons can be taken in account by a displacement 
of the interface to which the new reflection coefficient 
r' is referred with respect to the upper crystal plane to 
which reflection coefficient r refers in formula 01 After it 
we can rewrite 0] in the form 

Zb -1 ^ tan(fcoa/2) 



r exp(— j^) 



Zp 



(5) 



+ " tan{qa/2) 

The phase shift $ is rather small for qa < 1^^. 

If max{qa/2, koa/2) <C 1 the sine and tangent func- 
tions in ^ and [S] can be replaced by their arguments. 
Then we obtain from^Zs — ka/q that gives the approx- 
imation of the continuous medium. Really, if we com- 
pare Zb — ko/q = 1 / y^EcsHctf with Owe obtain — 1, 
q = ky/Ecs and El becomes the Fresnel formula for the 
refiection by a continuous dielectric half-space. The non- 
locality of the effective permittivity extracted through R 
in this case is negligible. The practical condition of ap- 
plicability of this approximation follows from numerical 
properties of tangent function and reads as 



a a , 
max I -, - — I < O.OI. 

A Aoff ^ 



(6) 
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This condition gives the upper frequency bound up to 
which formulas 121 and 13 are compatible with approxima- 
tion of the continuous medium. The frequency region El 
is the same for which transition layers in the theory^^ 
have no visible influence. For practical implementations 
of MTM the frequencies satisfying are not interesting. 

At frequencies ^ the normalized impedance Z = 
•v/scff/Mcff can be understood through parameters of a 
periodically loaded transmission line (TL). Formula [51 
obtained as the strict solution of the boundary lattice 
problem is known in the theory of periodically loaded 
TL. It represents the so-called Block impedance of a line 
loaded by shunt lumped impedances with period a (for- 
mula (5.117) fromSi). The Bloch impedance can be de- 
fined as the characteristic impedance of the homogenized 
periodically loaded line^^ normalized to the impedance 
of the host matrix. Notice, that the dispersion equation 
of a TL with shunt loads is known (e.gi^): 

j 

cos{qa) — coa{ka) H '■ sin(fca), (7) 

2^1oad 

where k is the matrix wavenumber and Zioad is the 
normalized impedance of loads. The equation 1291 coin- 
cides with equation (5.226) fromSl, thoroughly derived 
for an infinite lattice of electric dipoles at frequencies 
ka < TT. So, the problem of the reflection and transmis- 
sion in dipole lattices can be correctly solved replacing 
the crystal planes by effective lumped loads and using 
the model of a periodically loaded TL. The reflection co- 
efficient i? of a semi-infinite lattice is directly related to 
its Bloch impedance Z = Zb- If we measure R and T 
for a composite slab with integer number of cells across 
it, these parameters will uniquely determine q and Z 
(and, consequently, non-local MP). Thus, one can con- 
clude that MP calculated through R and T using (21 and 
l^lare transmission-line material parameters (TLMP), in- 
troduced for TL networks in 2002 in worksiSiSi. The idea 
of TLMP can be better understood from the comparison 
of Figs, n and 12 

To understand the reason why TLMP are non-local let 
us show that Zb can be expressed in terms of the Bloch 
expansion for microscopic fields in lattices. Namely, 
if a linearly polarized eigenmode with wave vector q 
propagates along the axis z in a lattice with periods 
{bx, by, bz = a) the Bloch impedance is equal to the ratio 
of zero-order spatial harmonics denoted as Eq and Ho in 
expansions: 

where Gn = (2nnx/bx, 2'Kny/by, ) are multiples of 
generic lattice vector. FromjHlwe have 



where the index TA means the transverse averaging that 
allows one to get rid of Gx, Gy in (HI It is evident 
that these integrals are exactly equal to Eta{—o-/2) — 
ETA{a/'2)exp{jqa), HTA{-a/2) = HTA{a/2) exp{jqa). 
The ratio Zb ~ i?TA(±a/2)/77i?TA(±a/2) can be inter- 
preted as that of effective voltage and effective current 
taken at the input or output of the cell. It is an al- 
ternative definition of the normalized Bloch impedance 
of periodically loaded TL with unit cell length oi. The 
Bloch impedance that determines (together with q) non- 
local MP of a lattice is equal to: (a) the characteristic 
impedance of the equivalent homogenized TL, (b) the ra- 
tio of fundamental Bloch harmonics of electric and mag- 
netic fields. The last fact is the reason of the non-locality 
of TLMP. 

Really, in the definition of local MP we deal with av- 
eraged electric and magnetic fields < E >, < H > 
which are simple integrals of the microscopic field over 
the domain around the observation point. Integrat- 
ing Maxwell's equations for microscopic fields does not 
change these equations and they hold for < E > and 
< H > and averaged polarizations. However, they do 
not hold for zeroth Bloch harmonics of microscopic fields 
taken separately. All high-order terms of expansionj^lgive 
a nonzero contribution into < E > and < H >! Treat- 
ing Eo and Hq substituted into Maxwell's equations as 
averaged fields leads to errors which can be interpreted 
as fictitious magnetization of p-lattices and fictitious po- 
larization of m-lattices. These fictitious lattice responses 
have no physical meaning, consequently, one should not 
expect the locality of TLMP. 



III. EXTRACTION OF LOCAL MATERIAL 
PARAMETERS 

In order to explain the extraction of local MP from R 
and T of a MTM slab, let us, first, consider an infinite 
lattice of electric and magnetic dipole particles. Assume, 
that a plane wave with wave vector q = qza propagates in 
an infinite lattice shown in Fig. ^ The lattice is formed 
by orthorhombic unit cells of sizes bxbx a (for simplicity 
one can skip the case of the anisotropy in the (x — j/)plane 
which gives no effect for the 2-propagation). 

In the left panel of Fig. ^ the particles are assumed 
possessing both electric and magnetic moments. In the 
right panel the electric and magnetic scatterers are differ- 
ent particles. The formulas of the local field method are 
the same for both these geometries. This is so since the 
quasistatic interaction between electric (p) and magnetic 
(m) dipoles is absent in both these structures (any p- 
dipole is not affected by m-dipoles lying in the same crys- 
tal plane and vice versa). Let a reference (0- numbered) 
p-m-particle located at the origin have electric and mag- 
netic polarizabilities Oee, CLmm- These polarizabilities can 
be resonant in the same frequency range. Let £'1°'= and 
H^°^ denote the x-component of the local electric field 
and y-component of the local magnetic field, respectively. 
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FIG. 1: Presentation of a p-m-lattice as a set of dipole crystal 
planes. Left: every particle has both electric and magnetic 
moments. Right: electric and magnetic particles are different. 
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FIG. 2: Presentation of the p-m-lattice shown in Fig. |21 as 
(1) periodically loaded transmission line, (2) homogenized 
transmission line with Bloch impedance Zb and wavenum- 
ber q {k and rj are unloaded wavenumber and characteristic 
impedance, respectively), and (3) effective magneto-dielectric 
medium with non-local MP Sef / and fi^f / . 



Then the p-dipole and m-dipole of the particle can be ex- 
pressed as: 



Po = aeeE^°'^, niQ = amniH 



loc 



(10) 



Local fields are produced by other particles located at 
points {x — Tixb, y = Uyb, z = riza)) and can be expressed 
through Po ^iid mg since p(n) = pQexp{—jnzqa) and 
TO(n) = moexp{—jnzqa): 

j^loc ^ CeePO + Cemmo, = Cmmmo + CmePO- (H) 

Substituting ITTl into [Till we come to a system 



Po = CeePa + Ce,„mo, 



1 



■Too — Cmm'niQ + CmePOi 



(12) 
(13) 



A closed-form expression for cross-coupled interaction 
factor Cem = Cme was derived ir^. For Cee = Cmm/'rf 
one can find in^'' the expression: 

rjLO f sinfca \ . k'^ 

2,0^ \ cos ka — cosqa J oirem 



(14) 



where k = uj^^mSm is the host matrix wavenumber and 

V = \J P-ml^m is its wave impedance (here Em and /i„i are 
absolute permittivity and permeability of the matrix). In 
I14l the trigonometric term describes the wave interaction 
of the reference p-dipole with p-dipoles of other crystal 
planes, the real value C^p describes the near-field inter- 
action with p-dipoles of other crystal planes. The remain 
of El is the real value Co which is responsible for the 
near-field interaction of the reference p-dipole with other 
p-dipoles of the reference crystal plane. For Co a simple 
expression was derived in-^: 



1 / cos khs . , , 



where s « 0.6954. At very low frequencies kh < 0.5 the 
approximation Co ~ 0.36 is quite accurate^^. 

Equating the determinant of El El to zero and ne- 
glecting Cnf (as it was shown in^^, it is a suitable ap- 
proximation for the case a < 6) we obtain the dispersion 
equation (an alternative equivalent form of it was derived 



cos qa = cos ka — sin ka 



G 



X 

T 



C X 



I SIT? ka\ 



GX 



4 4 y 4 

Here the notations were introduced: 



■ sm qa. 



koa 



e^nV 



(15) 



(16) 



Co 



In El y = ob^ is the unit cell volume. In the lossless 
structure the fundamental relations hold for electric and 
magnetic polarizabilitie a^i^i^ : 



Im- 



1 



fc3 
67re„ 



Im- 



1 



67r/i„ 



(17) 



Then the parameters G and X in El turn out to be real- 
valued in lossless lattices. 

Consider the averaged fields < E > and < H > defined 

as 



<E> {r) = ^ j E^{v + r')dPv' 



<H>{r)^^ j Hy{r + v')d 



(18) 



(19) 



Here r' lies within the unit cell and E^, Hy are com- 
ponents of the microscopic fields. For microscopic fields 
Maxwell's equations hold: 

rotE = -juj{fimtl + M), rotH = jw(£,„E + P). 
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Integrating these equations around the observation point 
with the use of |S1 we come to Maxwell's equations for 
averaged fields and polarizations: 



q < H >^ ujEm < E > +UJ < P >, 



(20) 



q < E ujfim < H > +LJ < M >, (21) 
whose solutions take form 



<E>, <H > 
<P>, <M > 



El, Hl 



Pl 



Mr 



cxp(-jgz). (22) 



The amplitudes El, Hl are related with Bloch harmon- 
ics En, Hn of transversely averaged fields as 



El 
Hl 



Eo/qa 
Ho/qa 



E 

n=l 



En/{qa - 
Hn{qa 



2nn) 
2Tm) 



(23) 



and at the reference particle center a; = y = z = Owe 
have < P >= po/V and < M nio/V. Equations EO) 
and 12 II at z — are equivalent to relations 

2 2 7 _ < -g > (0) 

, =.eL,.L, Zl = ^j^^=^-. (24) 

Here coefficients el and fj.L are defined by formulas 
<P>(0) f ieL-eo)<E>{0) 

(25) 



<M> (0) i (AiL - Mo) < > (0) 



and have physical meaning of absolute permittivity and 
permeability. They should be local even at moderately 
low frequenciesHsince relate averaged fields and averaged 
polarizations instead of zero-order Bloch harmonics. It is 
difficult to prove strictly this fact (see also the discussion 
inffi), but the explicit examples confirm it. 

Let us assume that we know parameters G and X de- 
scribing the crystal plane response. The evident way 
to find local MP is to express polarizabilities age, Qmm 
through G, X using formulas 1161 and then to apply the 
well-known Clausius-Mossotti formula expressing El via 
ttee and ^L via amm- However, this formula is quasi-static 
and is hardly accurate at moderate frequencies. Formu- 
las 12] show another (not static) way for obtaining local 
MP through G and X. Really, from [El [H we have 



sin ko 



G cos qa — cos ka 



2 

TOO I — 



sin ka 



X cos qa — cos ka 



= Too 



Wo 



sm qa 



cos qa — cos ka ' 



smqa 
cos qa — cos ka 



These relations give the auxiliary coefficient one needs in 
order to find Zl- 



Too 



?? < P > (0) 
< M > (0) 



' G / cos qa — cos ka 



G sin ka 



X y cos qa — cos ka — 



X sin ka 



(26) 



From 123 we obtain (see also 



^29.30. 



7 = vZl 



(gL - £o) 

(ml - mo) ' 



However, from 1241 one has 



£l = 



Ml = —Zl- 



ujZl 



(27) 



Therefore the wave impedance Zl can be expressed as 



^k - 



(28) 



So, from G, X one can find q, then 7, then Zl, and 
finally local MP. 

Practically, what is suggested is the generalization of 
the quasistatic approach based on the Clausius-Mossotti 
formula to the range of moderate frequencies. In the 
quasistatic limit our result for el and ml must be the 
same as this classical result. Let us prove it for simplicity 
for a cubic p-lattice. In this case X = Q and [151 simplifies 
to 



cos{qa) = cos(fca) 



G 



sin(fca) 



(29) 



Substitution of 7 — > 00 into 1281 gives Zl = rjk/q and 
q — fco\/£L- For low frequencies {ka <C tt) we take in 
account two terms of the Taylor expansion of trigono- 
metric functions in 1291 and substituting expression 1161 for 
G obtain: 



1- 



{koa)^£L 
2eo 



(fcoa)^£„ 
2eo 



{koa)^£rn ( 



1 



(fcoa) e„ 
6eo 



ye„,Re 



Co 



Neglecting the small term (fca)^/6 <C 1 we obtain the 
Clausius-Mossotti equation for local permittivity: 



£L 



1 



(i)-& 



VRe 



Our factor Co = 0.36 differs slightly from classical 0.33. 
This difference is related with the ignorance of near-field 
interaction between adjacent crystal planes. The needed 
correction can be easily introduced in the theory. Taking 
in account three terms in the Taylor expansion of cosine 
and sine functions gives the frequency corrections of the 
order {koa)"^ and {koa)^ to the Clausius-Mossotti rela- 
tion. These corrections were obtained in another way in 
work--. A similar proof can be done for p — m— lattices 
where one comes in the limit case (ka, qa) <C tt to the 
Clausius-Mossotti formulas for both el and ml- 
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Now, let us relate e^, with R and T coefficients of 
a finite slab. The simplest way is to find q and Z — Zg 
entering 13 inverting El in a usual way and then find G 
and X through them. With G, X we can calculate local 
MP. Let us then relate Bloch impedance Zb with G and 
X. For it one can apply the transfer matrix method. 
For the case under consideration (the normal propagation 
with respect to (x — y) crystal planes) the transfer matrix 
of a lattice unit cell F with components {A, B, C, D) is 
defined as 



Eta (- 
Hta 



A B 
G D 



Eta (f ) h 
. Hta (!) . 



(30) 



It can be calculated as a product of transfer matrices of 
its three portions: two a/2-long pieces of the TL mod- 
eling the host medium and the lumped circuit modeling 
the p-m-crystal plane (see Fig. EJ: 



Frf,F 



(31) 



Here 



■ TL 



Values G 





'1 Z' 


, Fyz = 


Y -1 



■cos(f) isin(^; 
[jsin(f) -cos(%)J 

(32) 

-jY and X = —jZ are parameters we in- 
troduced above as G and X and are given by^J For a 
lossless structure G and X are equal to the shunt imit- 
tance and series reactance of the loading, respectively. 
Parameters Y and Z are known in the theory of pla- 
nar grids as the grid electric admittance and the grid 
magnetic impedance and describe the ratio between tan- 
gential electric and magnetic fields in the grid plane and 
electric and magnetic surface currents induced in the grid 
(all averaged over the grid periods). Even if particles 
have sizes comparable to the period a one can attribute 
effective currents to the crystal plane z — 0. This proce- 
dure is basically the same as replacing the finite particles 
by point p- and m-dipoles which is an excellent approxi- 
mation for many practical cases. 

The Bloch impedance is given by relationSS: 



B 



sm qa 

After some algebra we obtain from 13 11^ 



Zb 



sm qa 



Gcos^ {ka/2)- Xsm%ka/2) + sin ka 



(33) 



(34) 



If one knows q and Z one can express G and X through 
them using equations 1151 and 1341 Then one finds local 
MP through G and X using [771 and 

To conclude this section one can notice that studying 
the transfer matrix of a cascade of lattice unit cells it is 
possible to prove that TLMBi^i^iSi are not mesoscopic. 
Being found from R and T of the layer containing N lat- 
tice unit cells across it these MP can be applied for slabs 
with arbitrary thickness (if divisible by a). Therefore, 
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(b) 



FIG. 3: The lattice of pairs of n-particles |(a)| and the lattice 
of in frared SRRs suggested by J.B. Pendry and S. O'Brien in^ 

m 



one can attribute TLMP found for finite slabs to infinite 
lattices. Considering the unit cell transfer matrix F one 
can also prove that constitutive parameters suggested in 
worksi2*i& and discussed ini& are equivalent to TLMP. 



IV. RESULTS 

To illustrate the locality of one set of MP and non- 
locality of the other one, two numerical examples of the 
direct calculation of them through the known polariz- 
abilities of particles are presented. The algorithm of this 
calculation is as follows: 



(«. 



ee 1 ^mm 



)~-{G,X) 




{EcS, Mcff) 
{£L,fiL) 



(35) 



The first arrow implies formulas 1 161 the second one im- 
plies the solution of 1 151 and the use of formula IS^ in the 
upper case and|^in the lower case. To find local MP for- 
mulas |5BH2H1 were used. The first numerical example cor- 
responds to a cubic lattice formed by capacitively loaded 
wire dipoles. The second example corresponds to a lat- 
tice of pairs of fJ-particles as is shown in Fig. |3| (a). In 
the first example the polarizability of particles is nearly 
static flee ~ I'^Gq, where I is the effective length of a wire 
dipole and Co is the loading capacitance. However, with 
rather large Co the effect of the presence of non-resonant 
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FIG. 4: |(a)| -dispersion in a lattice of non-resonant dipoles; 
|(b)| -dispersi on in a lattice of resonant p-m-particles (pairs of 
sTparticles) ; |(c)| -non-local MP for the lattice of non-resonant 
dipoles; |(d)| -non-local MP for p-m-particles; [(f)] -local MP for 
the lattice of non-resonant dipoles; |(f)| -local MP for the lattice 
of resonant p-m-particles. 
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FIG. 5: |(a)| -refr action index extracted in— from R, T coeffi- 
cients of a slab comprising the lattice of silver SRRs with pe- 
riod a — 600 nm; |(b)| -normalized Bloch impedance extr acted 
in^°; |(c)| -non-local permittivity extracted for this lattice; [(d)] - 
non-local permeability; |(e)| -local permittivity; |(f)| -local per- 
meability. 



dipoles can be significant. This effect is a wide stop-band 
with staggered mode near the first lattice resonance as 
one can see in Fig. 0] (a) . Notice, that quahtatively same 
results as shown in Fig. ^ (a,c,e) should correspond also 
to lattices of metal spheres of radius 2 — 4 mm which 
also behave below 15 GHz (the lattice period was cho- 
sen a — \Q mm in both examples and the host matrix 
was free space) as non-resonant dipoles with rather high 
static polarizability. A lattice of electric dipoles cannot 
have local magnetic susceptibility. And the relative local 
permeability ^rL is identically unity in Fig. 01(e). The 
locality conditions are satisfied until the 1st lattice reso- 
nance which occurs at 7.5 GHz. Non-local MP in Fig. 01 
(c) contain non-trivial permeability and imaginary parts 



of £off and /icff have opposite signs. Even below 7.5 GHz 
we observe in Fig. 01 (c) the violation of causality. 

In the second example both electric and magnetic po- 
larizabilities of doubled fi-particles are resonant. We cal- 
culate them using formulas22i^: 



A 



jujT 



B 



(36) 
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where Lq is the inductance of metal rings, / is the effec- 
tive length of the electric dipole induced in the particle, 
R is the ring radius shown in Fig. O In this example the 
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p-m-particles resonate at 6 GHz. Loss factor T was taken 
negligibly small in order to avoid the mess in complex so- 
lutions of equationEl The second example is illustrated 
by plots in Fig. El(b,d,f). The staggered mode holds in 
the narrow band 5. 8... 6 GHz. In this band local MP have 
no physical meaning. Only within this band the locality 
of ErL, ^J■rL IS lost. Outsidc it (cvcu at higher frequen- 
cies) the sign of Im(erL, Mi-l) is correct and the frequency 
behavior of Re(erL,MrL) in the low-loss region is causal. 

In our third example we study the structure for which 
the input data are R and T coefficients of a slab. It is 
filled by a lattice of silver split-ring resonators (SRRs) 
shown in Fig. |31 (b). This structure operating in the 
infrared range was studied ir&iS. For this case the algo- 
rithm is as follows: 

{R, T) ^ (q, Zb) ^ (G, X) ^ /^l). (38) 

In algorithm 1381 the first arrow corresponds to the inver- 
sion of formulas|21and|21 However, q and Zb were already 
found in^-^" and one can use the data presented in P. 56 
o£^ for n = q/k^ and Zb- In the second step of 1381 one 
can use formulas 1151 and 1341 solving them as an equation 
for G and X and finally find local MP from 1^1^ 

To calculate the MP in the third example the geome- 
try of infrared SRRs are not important since the input 
data are values of R and T. The period a in this case is 
equal a = 600 nm and the time dependence is exp(— iwf), 
as it is adopted in optics. The input data taken fromifl 



are shown in Fig. [S] (a,b) and the results are presented in 
Fig. |Sl(c-f). The difference between the results obtained 
for local and non-local material parameters confirms the 
theory developed in the present paper. Parameters ex- 
tracted by a direct inversion of Fresnel formulas [5] and 
|31 are non-local and parameters obtained using the sug- 
gested theory satisfy to the locality conditions. 



V. CONCLUSION 

In this paper it is demonstrated that local material pa- 
rameters of lattices can be introduced not only at very 
low frequencies but also in the region of moderately low 
frequencies where MTM operate. These local MP com- 
plement a pair of non-local MP which were obtained be- 
fore in known works. Though at moderate frequencies 
local MP allow to solve boundary problems for MTM 
slabs only at expense of introducing two transition lay- 
ers at the sides of the slab, these MP are very important. 
Only these MP do not depend on the incidence angle and 
only they can be used to study the interaction of MTM 
with wave packages and evanescent waves. It is then im- 
portant to learn to extract these parameters from R and 
T coefficients of MTM slabs. The simplest algorithm of 
this extraction is explained in the present paper, and the 
numerical examples show that the obtained MP satisfy 
to locality conditions. 
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